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Abstract 

Interfacial profiles and interfacial tensions of phase-separated binary mixtures of Bose-Einstein 
condensates are studied theoretically. The two condensates are characterized by their respective 
healing lengths £i and £2 and by the inter-species repulsive interaction K. An exact solution to 
the Gross-Pitaevskii (GP) equations is obtained for the special case £2/^1 = 1/2 and K = 3 / 2 . 
Furthermore, applying a double-parabola approximation (DPA) to the energy density featured in 
GP theory allows us to define a DPA model, which is much simpler to handle than GP theory 
but nevertheless still captures the main physics. In particular, a compact analytic expression for 
the interfacial tension is derived that is useful for all £1,^2 and K. An application to wetting 
phenomena is presented for condensates adsorbed at an optical wall. The wetting phase boundary 
obtained within the DPA model nearly coincides with the exact one in GP theory. 
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I. INTRODUCTION 


Phase separation in gas mixtures goes against our intuitive notion of the entropy of 
mixing. However, at ultralow temperature where quantum physics rules, this possibility has 
been realized experimentally, in particular in binary mixtures of Bose-Einstein condensates 


(BECs). See, e.g., [l| for an overview. Let us recall briefly the state-of-the-art research in 
this field, paying special attention to experiments. 

At the beginning of this century weakly phase-segregated binary Bose-Einstein systems 

71 More recently, strong phase separation has been 
12 ]. Moreover, also in an ultracold dual component 


have been observed experimentally Of 
realized by various research groups B- 


gas of thermal atoms spatial separation of components has been achieved 


13), while many 


more degenerate Bose mixtures were produced in which phase separation is possible 14B16], 


The theoretical and experimental physics of multi-component condensates is well explained 
in Refs, [h], 18]. The statics and dynamics of phase separated BECs have been extensively 

The interfacial phenomomenology in Bose mixtures has been 


addressed in Refs, 
explored in Refs. 
Refs. 


Hail. 


32 


40 


39) and the phase diagram at finite temperature was investigated in 


42], 


Our focus in this paper is on the calculation of static interfacial properties of BEC binary 
mixtures within Gross-Pitaevskii (GP) theory 43). Research on this problem, in particular 
on the interfacial tension of phase-separated mixtures of BECs 24), has led to interesting 


results 


23 


44). In particular, in [44] accurate analytical approximations (e.g., series ex¬ 


pansions) have been provided, covering certain ranges of condensate healing lengths and 
interparticle repulsive interaction strengths. While these results are useful, there is still a 
need for i) more exact solutions within GP theory, and ii) a simpler model which can pro¬ 
vide a compact and insightful expression for the interfacial tension that can be used over 
the entire parameter range. Our aim is to contribute advances to meet both of these needs. 

This paper is organized as follows. Section II deals with the frame-work of GP theory. 
Section III is concerned with exact solutions to the GP equations describing interfacial 
profiles and announces our exact solution for a special choice of physical parameters. Section 
IV defines our simpler model through the so-called double-parabola approximation (DPA) to 
the GP Lagrangian and presents the solutions for the interfacial profiles within this model. 
Section V treats the application of the DPA to the interfacial tension and presents our 
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compact analytical expression for this quantity. Section VI is concerned with the application 
of the DPA to the wetting phase transition that can occur when the condensates are adsorbed 
at an optical wall. This section also clarifies how the exact solution announced in section II 
could be found by virtue of insights gained through the DPA. Finally, our conclusions and 
outlook are given in Section VII. 

II. GROSS-PITAEVSKII THEORETICAL FRAME-WORK 

A. Gross-Pitaevskii Lagrangian 


We start from the Lagrangian L and action S' of a two-component BEC 

S'[T 1 ,'F 2 ] = J d tL = I dtdr£ 

with the mean-field Gross-Pitaevskii Lagrangian density 45|, |46] 


ih 


C(*1, *2) = - S(T 1; T 2 ), 


( 1 ) 


( 2 ) 


3 =1 


in which the Hamiltonian density is 


2 r 


S(T 1 ,T 2 ) = 

3 = 1 


A_|v>P J f + M|«P j | 1 + 3 l 2 |l' 1 | 2 |<I < 2 | 2 


( 3 ) 


where, for species j , Tj = \l/j(r, t) is the wave function of the condensate playing the role 
of order parameter; rrij is the atomic mass, g 31 = 47t h 2 ajj/m,j > 0 is the strength of the 
repulsive intra-species interaction, g\ 2 = 27r/i 2 ai 2 (l/mi + l/m 2 ) > 0 is the strength of the 
repulsive inter-species interaction and cijj> is the s-wave scattering length, relevant at low 
energies. 

By introducing the following dimensionless quantities, s j = rwith ^ = h/ yj2mjrijogjj 
the healing length and rijo the number density of condensate j in bulk, Tj = t/tj, Vb = 
Tj/y'njo, and K = gi 2 /y/gng 22 , where tj = h/fXj, and fij = gjjfijQ the chemical potential of 

condensate j, we scale the Lagrangian density in (J2J) and Hamiltonian density in (j3J) to 

2 

i 1 ~ 

( 4 ) 


C 


£(^1,^2) = 77^ = V r -£(^1,^2), 

2P 0 J=] 2 


with 


£(^ 1 ,^ 2 ) = 


£ 

2D 


E 

3 =1 


|v s v/ + % 


14 1 


+ K\^\H 2 , 


( 5 ) 
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where the pressure Pq is given by fi 2 / 2g 3J , which takes one and the same value in both 
condensates at two-phase coexistence. Next we make a transformation of the dimensionless 
Lagrangian density by writing 


lM s . pTj) = 0j( s ,> r j) e " 


( 6 ) 


We then have a Lagrangian density in terms of the new order parameters 0,-, 


2 r • 

i 


£ (0i, fa) = £ {fae lT \fae ZT2 ) = b - fad Tj faj) ~ | V s .0_ 


3 = 1 


in which the potential V takes the form 


V(0i,0 2 ), (7) 


V(01,0 2 ) = X] 

J=1 


H0j| 2 + 


10 / 


+ K\fa\ 2 \fa' 2 


( 8 ) 


We recall that when K > 1, the two components are immiscible and a two-phase segregated 
BEC is formed 24]. 


B. GP Equations 

By considering a variation 0* —i► fa- + <50* and requiring 5S/S(p* = 0, we obtain the 
Euler-Lagrange equations 

<9£ <9£ . . . 

o i * /,■ 7)77) Tn’ 7 £2, (9) 

<90* 3 <9(9 r .0*) 

yielding the time-dependent GP equations 


72 , . . 


id Tj4>j = -V s .0,- + —; j = 1, 2 


( 10 ) 


or 


= [-V^ - 1 + 10/ + A-I^rj j = 1,2 (j ^ /). (11) 

Note that these reduce to the time-independent GP equations (TIGPE) when the order 
parameter (pj(sj,Tj ) = 0jo( s j) is time-independent (so that fa(sj, r^) = 0jo(Sj)e _ * Tj ' is sta¬ 
tionary), namely 


<9V 

V s>tO = 7T7J - j J = £2, 


90 


( 12 ) 


jo 


which leads to 


— Vg — 1 + 10jo| 2 + K 10j v o| 2 0jo — 0; j — 1, 2 (j ^ j') 


(13) 
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C. Boundary conditions for interface profiles 


For describing a static planar interface at z = 0, separating bulk condensate 1 residing 
in the half-space z > 0 and bulk condensate 2 residing in the half-space z < 0, we limit our 
attention to order parameters that are translationally invariant in the x and y directions. To 
keep the notation simple, we will not change the name of the function: (frjo(sj) —y 4>jo{Pj), 
where pj = z/£j. For describing an interface the TIGPE must be solved with the boundary 
conditions 


0io(pi -t oo) = 02 o(P 2 ->■ -oo) = 1 
02 o(d 2 -)• oo) = 010 (pi -oo) = 0 


(14) 


III. EXACT SOLUTIONS FOR INTERFACIAL PROFILES 


A. The strong segregation limit K —> oo 


In the interesting limit of strong inter-species repulsion K —> oo, the segregation (i.e., mu¬ 
tual exclusion in space) of the phases becomes complete. Numerical solution of the TIGPE 
indicates that the overlap of two order parameters becomes zero and that the interaction 
term 7U| 0 X | 2 10 2 1 2 in the potential (JSD becomes negligible. Consequently, the GP equations 
decouple in this limit (in spite of the divergence of K). In this limit the simple exact solution 
to the GP equations for the interface consists of a pair of adjacent “tanh” profiles 


0jo(/b) = tanh 


(- 1 ) 


3 +1 Pi 

y/2. 


; 3 = 1,2. 


(15) 


The interface position is conveniently marked by the common point of vanishing of the 
two order parameters. Note that the interface consists simply of two adjacent “hard wall” 
profiles [47|. It is noteworthy that a (single) condensate wave function at a hard wall is 
mathematically similar to the Ginzburg-Landau (GL) order parameter for a strongly type-I 
superconductor at a normal/superconducting interface in the limit k —> 0, with k the GL 


parameter 


48]. 


At this point we make two remarks. Firstly, the fact that the interaction term iF|0i| 2 |02| 2 
in the potential (J8]) becomes negligible, can be understood analytically using the DPA, which 
will be introduced in the next section. Secondly, the appearance of the “tanh” function is 
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not unique to the strong segregation limit. In the next subsections we will see two more 
examples of the occurrence of this function in exact solutions for finite K. 


B. The exact solution of Malomed et al. 


Interestingly, for K = 3 and for the symmetric choice £1 = £ 2 , an exact solution to 
coupled differential equations mathematically identical to the GP equations was provided 


by Malomed et al. 


53| in the physical context of domain boundaries in convection patterns. 


The solutions are surprisingly simple and again involve the “tanh” function, 


no — 


1 + tanh 


\/ 2 £ 


; <j>20 — 


1 — tanh 




(16) 


We will come back to this special case when we discuss the interfacial tension in section 
V. There it will become clear that the DPA can explain why we might have anticipated that 
the exact solution is simple for this special case K = 3 and £1 = £ 2 . 


C. An exact solution for an asymmetric case (£1 / £ 2 ) 


We have uncovered another exact solution. It applies to the case £ 2 /£i = 1/2 and K = 3/2 
(the roles of £1 and £2 can be interchanged). The solutions to the GP equations are now the 
following profiles, with £1 = £, and £2 = £/2, once more involving the “tanh” function, 


010 



020 


1 

2 


1 — tanh 



(17) 


The heuristic procedure through which this serendipitous solution could be found has been 
strongly guided by insights provided by the DPA. We will discuss this in detail in section VI 
when we treat the wetting problem within DPA. In this and the preceding two subsections 
we have disclosed that the DPA is useful for guessing exact results. We now turn to the 
precise definition of the DPA in the context of GP theory and to the study of the remarkable 
properties of the simple model defined through the DPA. 
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IV. DOUBLE-PARABOLA APPROXIMATION (DPA) FOR INTERFACIAL PRO¬ 
FILES 


The idea of approximating a double-well potential by a piece-wise parabolic function is 
old and generally dictated by the desire to work with piece-wise harmonic potentials that 
allow one to solve the equations of motion exactly using simple functions, the behavior of 
which is easy to interpret physically. In the context of surface and interfacial phenomena one 
of the first to implement this approximation for the bulk free-energy density was Hauge 49]. 
We follow this line of thought and expand the quartic potential V in (jHJ) about its (local) 
minima, which correspond to the bulk values for the order parameters. For obtaining the 
interface profile in the half-space z < 0 we make use of the expansion about bulk condensate 
2 and for the half-space z>0we expand about condensate 1, in the following manner, 


\<j)j\ = 1 + 6j] \(j)jt\ = Sji, with 


*>0, (j,j r ) — (1, 2) 


(18) 


^< 0 , = ( 2 , 1 ) 

in which the real numbers 6j and 5j are treated as small perturbations. We expand the 
potential up to second order in ej and 5j and so arrive at two “quadratic” potentials, each 
of which is to be used in the appropriate half-space, 


Vdpa(0i ,<h) = 2 (\H - l) 2 + (K - 1)1 


b>o, (hi') 

[z<o, (j,f) 


( 1 . 2 ) 

( 2 , 1 ) 


(19) 


This defines the DPA for the potential energy density and can be interpreted as the potential 
for a model that is related to, but different from, the original GP theory and which we will 
call the DPA model. The following technical remark is in order: In view of the structure of 
the TIGPE we anticipate that, at the interface position z = 0, it will (in general) be possible 
to preserve continuity of the order parameter functions and their first derivatives. Indeed, in 
view of (fT2]l . continuity of the potential V but discontinuity of one of its (partial) derivatives 
with respect to the order parameter(s), at z = 0, will induce a discontinuity in the second 
derivative of the interface profile functions (ftjo(Pj)] j — 1, 2. This is a mild singularity, often 
imperceptible to the eye in an interface plot. Experience has taught us that this singularity 
has little or no effect on the qualitative features of the phenomena under study, provided 
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one stays away from bulk criticality or similar conditions, at which the (local) minima of 
the potential may merge or undergo some other drastic change. 


A. DPA for the GP equations 


The double-parabola-approximated GP equations are obtained by replacing the potential 
V in the Lagrangian density (J7]) by the Vdpa hr (fT9l) and deriving the Euler-Lagrange equation 
(12]), which leads to 


+ (K - l)4>? 


with 


2>0, (j,j') = ( 1,2) 

^<0, (j,j') — (2,1) 


(DPA), (20) 


where (f>j = \(j)j\ exp (idj). Note that, alternatively, these linearized equations can be obtained 
directly from the GP equations of (fTTT) by expanding the order parameters to first order in the 
perturbations about their bulk values, similarly to what was done in 24], For studying the 
interface structure, it would be sufficient to limit our attention to real order parameters and 
real perturbations. Nevertheless, we insist on deriving the DPA in complex function space 
because we have the intention to apply our results in future work to dynamical properties 
such as dispersion of phonon excitations and capillary wave excitations using the Bogoliubov- 
de Gennes formalism. 

When the order parameters are stationary, namely (j)j(sj,Tj) = (j)jo(sj), we obtain the 
DPA to the TIGPE, 


-VsA-o + 2 


jo 


5 )=o 


with 


+ {K ~ l)0j'o — 0 


^>0, (j,f) — (1, 2) 

^<o, a; 1 ) = (2,1) 


(DPA), 


( 21 ) 


where o r = |^ 0 | exp 

Note that the DPA equations appear decoupled in the order parameters. However, the 
form of the equations depends on which bulk phase is chosen as starting point for the 
expansion. Since the bulk boundary conditions are different on either side of the interface, 
the order parameters are implicitly coupled. Of course, the apparent decoupling greatly 
facilitates analytical calculations. 



B. DPA solutions 


For obtaining static interface profiles we can take all functions to be real, and independent 
of the coordinates x and y (translational invariance in the directions parallel to the interface) 
and independent of time. The solutions 4>j(pj,Tj) = <f>jo(pj) that solve the TIGPE with the 
suitable boundary conditions (fT4l) are simple exponentials. A unique interface is obtained 
by matching the solutions for z > 0 to the ones for z < 0 at z = 0 with the requirement 
that the functions and their first derivatives be continuous at z = 0. This leads to 


0jo(Pj) — 1 

fij'oiPj') = 


P 


ol T j3 




a 


ol + (3 


e~P\ p i’\ 


with 


2>0, (j,/) = (l,2) 

*<0, (j,f) — (2,1) 


(DPA), 


( 22 ) 


where a = \[2 and f3 = \/K —1. Note that Ao and Chui 24] already discussed these 
functions in the frame-work of perturbation expansions and identified //J as the penetration 
depth of condensate j into condensate j' for j ^ f. 

In the following two hgures we compare order parameter profiles calculated within the 
DPA model with numerically exact order parameter profiles solving the GPE (TT3h . For the 
symmetric case £i = £2 = £, Fig.l shows the numerically exact order parameters together 
with their DPAs for K = 2 (moderately weak segregation). For an asymmetric case 2^ = £ 2 , 
Fig. 2 shows the numerically exact order parameters together with their DPAs, also for 
K = 2. Note how the DPA differs from the exact solution in featuring steeper profiles, 
corresponding to a smaller interface width. In spite of this quantitative difference the DPA 
appears to lead to a qualitatively correct interfacial structure. 


C. Validity conditions 


Ao and Chui 24j] pointed out that, in order for these solutions (j22li of the linearized 
GPE to be self-consistent within perturbation theory, a condition corresponding to (fairly) 
strong segregation must be fulfilled, which requires the penetration depth of condensate j 
into condensate i to be smaller than the healing length of condensate i, for i — 1,2. This 
condition is 

e j /v / 2(^-l)<6; (23) 
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K = 2 , 6/6 =1 



FIG. 1. (color online) Interface structure for weak segregation. Reduced order parameter profiles 
c/)jO, j = 1,2, are plotted versus z/6 for K = 2 and the symmetric case 6 = 6- The numerically 
exact profiles (black lines; GP) and the double-parabola approximations (grey lines; DPA) are 
shown. 


When this condition is not satisfied, as can typically happen for weak segregation (K > 1), 
a term quadratic in the penetrating condensate order parameter dominates a term of first 
order in the deviation of the penetrated condensate order parameter from its bulk value 24], 
Under such circumstances it is recommended to go beyond the DPA as far as the compu¬ 
tation of order parameter profiles is concerned. However, for energy and interfacial tension 
computations we find that including this nonlinear term brings only modest improvement 
over the interfacial tension calculated within our DPA model, when the result is compared 
with the (numerically) exact result in GP theory, even in the weak segregation regime. 

Our strategy in this and forthcoming works is, and will be, to consider the DPA as a 
model in its own right, based on and defined by the specific potential energy density (TTUT) . 
and to explore its predictions. As a first example of this strategy we will, in the next section, 
calculate the interfacial tension within the DPA model and compare it with the interfacial 
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K = 2 , 6/6 =2 



FIG. 2. (color online) Interface structure for weak segregation. Reduced order parameter profiles 
(f>jO, j = 1,2, are plotted versus z/G for K = 2 and the asymmetric case 26 = 6- The numerically 
exact profiles (black lines; GP) and the double-parabola approximations (grey lines; DPA) are 
shown. 

tension within GP theory. As a second example we will, in section VI, apply the DPA model 
for the calculation of the wetting transition in adsorbed binary BEC mixtures, and compare 
it with the wetting transition in GP theory. 


D. Strong segregation limit 


In the limit of strong segregation K —> oo, we have a/(a + j3) —> 0 and (3/(a + (3) —$■ 1. 
The order parameters become 


(fijo(pj) = 1 - e a|Pjl ; (pj'o(pj') = 0, with 


\z>0, 

< 0 , 


{j,f) 

U,f) 


( 1 , 2 ) 

( 2 , 1 ) 


(DPA) 


(24) 


and we notice that there is no overlap of the condensates (and also no gap in between them; 
they touch at * = 0). Complete segregation is illustrated in Fig.3 for the symmetric case 
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K = oo, 6/6 =1 



FIG. 3. (color online) Interface structure in the limit of strong segregation, K —>• oo. Reduced 
order parameter profiles (f)jo, j = 1, 2, are plotted versus z/6 for the symmetric case 6 = 6- The 
exact “tanh” profiles (black lines; GP) and the double-parabola approximations (grey lines; DPA) 
are shown. 

6 = 6 - 

When we relax the strong segregation slightly so that the condensates incur a small 
but nonzero overlap, numerical analysis of the full GPE reveals that the interaction term 
K 10 1O | 2 10201 2 in the potential (J8]) of the Lagrangian density f7]) is negligible compared to the 
other terms. Indeed, in spite of the fact that K is large, the density overlap is so small that 
the product concerned is small. Now, the DPA allows us to establish analytically in which 
manner the interaction term vanishes in the limit K —> oo. Inserting the DPA solutions 
for the order parameters (l22|i in the interaction term and calculating its value at £ = 0, we 
obtain 

a>i«(0)|V 2O (0)| 2 = k (6y) K 76 0 (DPA)l (25) 

which quantifies the vanishing of the interaction term for strong segregation. It is therefore 
clear that, in spite of the diverging coupling constant, the interaction term can be safely 
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ignored in calculations pertaining to totally segregated condensates. Consequently, the 
GPE decouple in this limit. Thus, as we announced in subsection III.A, the DPA solutions 
can provide us with an analytical result, (1251) . allowing us to gain physical insight into a 
interesting property of the GP theory. 


V. DPA MODEL APPLIED TO THE INTERFACIAL TENSION 


A. Grand Potential and Interfacial Tension from the GP Lagrangian 


We define a dynamical grand potential as the grand potential in equilibrium but with 
time-dependent order parameters, 


!2[<Ai, 02] = 2-P 0 / drVV(^i,&), 


(26) 


where the integration is over the volume V of the system and the grand potential density is 
defined as minus the Lagrangian density in (J7J), i.e., 


VVW’i.cfc) = -C(<h,<h) 


(27) 


Consider an infinite planar interface at z = 0 and assume translational invariance along x 
and y. Without loss of generality, we assume the order parameters (j)j 0 for stationary states 
to be real and arrive at the grand potential for the interface, 

r \ 2 

£ (Mi°) +^(^10)020) , (28) 

- 3 = 1 


Q 0 — Q[(f) w ,(j) 20 ] — 2 P 0 A / d z 


where A is the interfacial area. 

To calculate fl 0 , we first derive a “constant of the motion”. We multiply the TIGPE in 
(fl3|i by d z (j)j 0 and add up the two equations. We then integrate over z and find 

2 r „ i i 

- K^l o = 1/2, (29) 


Z! o) +^ 2 o-2^ 4 o 

3 =1 *" 


where the constant 1/2 is obtained by considering the (bulk) boundary conditions for 
For order parameter profiles that satisfy the TIGPE we may substitute (|29|) in (1281) . Then 
Oo reduces to 

~ ’ " ( 30 ) 


!2o = 4/2,,4 ^ / dz (0 f! o r ,\ 2 - P„V. 

3 =1 J 
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The interfacial tension is defined as the excess grand potential per unit area, 


712 = j — = 4P 0 [ dz ( d Pih °) 2 ’ ( 31 ) 

j= i ' 

which remains after the bulk grand potential = —PqV of a homogeneous phase has been 
subtracted. Note that at bulk two-phase coexistence the bulk grand potentials are the same 
for each condensate. Note that expression (T3T1) is valid for solutions of the TIGPE. If we 
wish to evaluate (nonequilibrium) interfacial tensions in GP theory for profiles that do not 
necessarily satisfy the TIGPE we must use (I25|) . 


B. Interfacial tension within the DPA model 

Expression ()3T1) for the equilibrium interfacial tension is independent of the form of the 
potential V. Therefore, we obtain the same expression if we start from the model defined 
by the DPA potential (TT9|) and use profiles that satisfy (I2TT) . Note that there is no simple 
relation between the value of the GP interfacial tension and that of the one defined within 
the DPA model. They are equilibrium interfacial tensions for two different models. Note 
that our approach is fundamentally different from a trial-function approach in GP theory. If 
we wish to consider the DPA solutions as trial functions wthin GP theory, we must use ()28li 
and may not use ()3B . Doing so would lead to an approximation that is far less useful than 
our DPA model because, for example, it would lead to an interfacial tension that diverges 
in the limit K j. 1, which is physically unacceptable. 

In view of these considerations, an analytic expression for the interfacial tension within 
the DPA model is obtained by evaluating (l3Tf in the DPA profiles (1221 . We find 

P2) 

This compact expression is insightful. It shows that, in the simplified model, the contribution 
from each condensate is proportional to its healing length. Furthermore, the expression 
interpolates between the strong and weak segregation limits by means of a function that 
depends only on K. How this DPA interfacial tension compares to the GP interfacial 
tension, is the question to which we now turn. 
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FIG. 4. (color online) The reduced interfacial tension 712 /T 0 ? f° r the symmetric case (£1 = £2 = £) 
versus the inverse interaction strength 1/A'. Shown are the numerically exact solution (GP, black 
line), the double-parabola approximation (DPA, grey line) and the scaled DPA (red line) obtained 
by multiplying the DPA result by 2/3 to match the DPA and GP values at 1/K = 0. 


C. Comparison with exact results in GP theory 


1. Strong Segregation 


In the strong segregation limit K —y 00 , 712 is the sum of the two wall tensions. It can be 
calculated directly by substituting the stationary solutions for the strong segregation 
limit into (131 D . so we have 


712 — 7 wi + 7 W 2 — 


4^2 


^ 0(61 + 6 ), for K -A 00 , 


(33) 


where 7 is the wall tension for the j-th component 47]. Note that the wall tension (for 
a single condensate) is mathematically similar to the tension of a normal/superconducting 


interface in the limit of strongly type-I superconductors 


48]. 
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Comparison with 7 i 2 >PA ^ given in (l32l) is immediate since 


(DPA) _ (DPA) . (DPA) o /on (C , C \ C TS , 

7i2 — 7 wi + Iwz ~ 2v 2P 0 (6 + £2)5 f° r K — y 00, 


( 34 ) 


where 7 * s the wall tension for j-th component. So, the ratio between the DPA and 
the GP result is 


(DPA) o 

—-= —, for K —> 00 . 

712 2 


For large K we can expand 


(35) 


7<? pa) = 2 V 2 p„ 1 1 -^ + o (411K, + &), 


(36) 


and observe that the interfacial tension approaches its limit with a sqnare-root singularity 
(with diverging slope) in the variable 1/A". We can compare this with the leading terms in 
the large-A" expansion of the GP interfacial tension, derived in 44], 


p I 9 ^66 1 

712 = Po I —:-2.056 


6 + 6 


+ ° (76)) (fi+&i 


(37) 


It is noteworthy that the A '~ 1//4 singularity as well as its amplitude 2.056... are universal 
in the sense that they are to some extent independent of the details of the theory. In 
particular, these features are common to the GP theory of BEC [43| and the earlier Ginzburg- 


Landau (GL) theory of superconductivity [48|. I 11 particular, the amplitude 2.056... was 


first obtained by Mishonov (1.03 in his units) using GL theory ^]. Subsequently this was 


511 ] and later applied to BEC [44| using GP theory. 


confirmed and elaborated 

The DPA result (1361) does not capture the A ' -1 / 4 singularity and displays a K ~+ 2 sin¬ 
gularity instead. The difference in the manner the slope diverges near 1/K = 0 can be 
appreciated in Fig.4, in which both the DPA and the GP interfacial tension are plotted as 
a function of 1/K for the symmetric case £1 = 6 - I 11 this figure we observe that the DPA 
curve is similar in shape to the GP interfacial tension. If we reduce the DPA by applying 
an overall prefactor of 2/3, we obtain the red curve, which follows the GP curve surprisingly 
well. Note that the reduced DPA intersects the GP curve in one internal point, to which we 
will return shortly, after discussing the other interesting limit, K 1. 
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2. Weak Segregation 


531 ] in a different physical context, who obtained 


In the weak segregation limit K j. 1, the two condensates tend to merge. The total 
density of the two condensates is nearly constant and displays a small depression around the 
interface. In this limit 712 was calculated by Barankov |52j, and before him by Malomed et 
al. c _ j 

712 ' Pci ifh - ~^— p ° «,+& ' 

When — £, it simplifies to 

712 = 2 \/K — 1P 0 £- 

Considering the DPA, when we take the limit K \. 1, we find that 7^ PA ^ 
proaches 

7 £ PA) = 2VK^1Po^i + £ 2 ). 

Comparing the DPA with the GP result for the symmetric case £j = £, we have 


111 


(DPA) 

712 

712 


= 2, for K i 1. 


(38) 

(39) 
ap- 

(40) 

(41) 


We conclude that the DPA model describes the interfacial tension in the weak segregation 
regime qualitatively correctly, since it shares the correct square-root singularity at K = 1 
with the GP expression. This can also be appreciated in Figure 4. 


3. Half-segregation and the exact solution of Malomed et al. 

We already noted, when we discussed Fig. 4, that the reduced DPA for the interfacial 
tension intersects the GP curve in one internal point. This happens at K = 3, for the 
symmetric case £1 = £ 2 - Interestingly, at this point in parameter space an exact solution 
to the GP equations was provided by Malomed et al. 53]] - We have already recalled this 
solution in (fl6|) . Note that the order parameters are perfectly symmetric. Not only are the 
healing lengths equal but the healing length also equals the penetration depth since \f2 = 
\JK — 1. The profiles cross precisely half-way their bulk values, at (f)jo = 1/2. Therefore, 
we denote this special case by “half-segregation”. 

It is instructive to observe that the DPA solutions (1221) are also perfectly symmetric in 
this case and display half-segregation. The interfacial tension within the DPA model is given 
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by 


712 


i? PA) = 2^2 >„£, 


while the GP value is precisely two thirds of this, 

4^2 


712 = 


77, 


(42) 


(43) 


VI. APPLICATION OF THE DPA MODEL TO WETTING PHENOMENA 


A. Wetting phase boundary 


In this section we apply the DPA model to the wetting phase transition predicted for BEC 


mixtures adsorbed at an optical wall, using the GP theory at T = 0 


M, 


55[. The wetting 


transition takes places when, e.g., a layer of condensate 2 intrudes between condensate 1 
and the optical wall (evanescent wave emanating from a prism). The condition for a wetting 
transition is the following surface energy equality, 


1W1 — lW2 + 7l2- 


(44) 


The curve in the (£ 2 /£i, l/AT)-plane that solves this equation is the so-called wetting phase 
boundary. In the hard wall limit (for a vanishing order parameter at the wall), the wetting 
phase boundary was established numerically |54] and an analytical solution was reported in 
the second paper of [44]. In this limit the wetting transition is of first order (discontinuity 
in the first derivative of the energy). The analytic solution for the phase boundary is 


T = 


72 


£1 _6 
L 6 7, 


(45) 


55|. It was 


For the more physical case of a softer wall, the wetting transition was studied in 
found that first-order wetting as well as critical wetting are possible. 

Our aim is to derive an approximate wetting phase boundary in the hard wall limit within 
the DPA model and to compare it with the GP result (l45lh To this end we first give the 
DPA for the wall tensions, 

7 ^ PA) = 272PoO, (46) 

as follows from (l32l) in the limit K —>■ 00 . Inserting these and our expression 
into (giD , leads to the wetting phase boundary within the DPA model, 


r (DPA) 

for 7)3 


\/K - 1 = -= 


72 L 6 


4_i 


. (DPA). 


(47) 
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FIG. 5. (color online) Phase boundary for the first-order wetting transition in adsorbed BEC 
mixtures, in the plane of inverse interaction strength and healing length ratio. The GP solution is 
shown (black curve) as well as the DPA (grey curve). 


Figure 5 shows the GP wetting phase boundary together with the DPA. Clearly, the two 
curves are almost coincident. Moreover, the DPA reproduces the parabolic character of the 
GP phase boundary near both endpoints, at 1/K = 0 and K = 1. 


B. DPA-assisted design of an exact solution to the GP equations 

In this subsection we explain how we found the exact solution (fTD) announced in section 
III.C and displayed in Fig. 6. We start by observing that the DPA intersects the GP curve 
precisely at £ 2 /£i — 1/2 and K = 3/2 (see Fig. 5). The asymptotic behavior of the order 
parameters can be read off from the DPA solutions, provided conditions (Pul) are satisfied. 
For z>0we can rely on the DPA since £ 2 /£i < a/ 2(A" — 1) = 1. An “up-down” symmetry 
occurs since £i/\/2 equals the penetration depth £ 2 /\/K — 1. For z < 0 we cannot rely 
on the DPA for the approach of 0 2 o towards 1, since £i/£ 2 > 1- In this case the approach 


19 









K = 1.5, 6/6 =0.5 



FIG. 6. (color online) Interface structure for an exactly solvable asymmetric case. Reduced order 
parameter profiles (f>jo, j = 1, 2, are plotted versus z/^\ for K = 3/2 and £i = 2£ 2 - The numerically 
exact profiles (black lines; GP), the exact solution (red dashed lines; analytical) and the double¬ 
parabola approximations (grey lines; DPA) are shown. 


towards the bulk density is governed by the decay length £ 1 /(2 \JK — 1), which is longer 
than the length £ 2 / predicted by the DPA. 

This information, together with the graphical observation that </> 10 takes a value of about 
l/\/2 at the point in space where </> 2 0 equals 1/2, suggest that the solutions to the GPE ought 
to be well approximated by the following skewed profiles presented in (fTTl) , which happen to 
solve the GPE exactly. Note that the ^-coordinate can be shifted so as to provide the profile 
crossing at z = 0, which facilitates comparison with the DPA which intrinsically features 
this position as the location of the interface. The required shift is 5 = arctanh(2 — y/E) in 
units of y/2^. Fig.6 shows these GP solutions together with their DPA counterparts. 

Incidentally, we note that the interfacial tension obtained for this exact solution 

2\/2 

712 = ^ (48) 
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is 2/3 of the value found within the DPA model for the same condensate parameters K = 3/2 
and £1 = 2£ 2 = f • 


VII. CONCLUSION 


In this work we accomplished two goals: i) We added an exact solution to the GP theory 
for interfaces in BEC binary mixtures. To our knowledge this is the first exact solution for 
an asymmetric system (with unequal healing lengths) and at finite inter-species repulsion 
strength K. We have been able to End this solution guided by information gathered by 
solving a simpler but related model, the so-called DPA. ii) We defined and developed the 
“DPA model”. We first derived the DPA for the potential energy density in the Lagrangian 
by expanding the order parameters about their bulk values and keeping the deviations into 
account to second order. Locating the interface center at z = 0, we next derived the DPA 
for the GP equations in each half-space, z < 0 and z > 0. The solutions and their first 
derivatives were then matched at z = 0. This led to unique simple analytical solutions that 
can be used efficiently to uncover and understand properties of GP theory. 

The power of the DPA model lies in its capacity to provide systematically analytical 
expressions for many physical quantities of two segregated BECs. We know of no other 
methods capable of doing this. Here, the excess surface energies at a hard wall and at the 
interface were evaluated within the DPA model. This provided a compact and useful expres¬ 
sion for the interfacial tension. As an application we derived the wetting phase boundary 
within DPA and obtained good agreement with the GP solution. Moreover, the DPA pro¬ 
vided crucial hints facilitating a successful guess of an exact solution to the GP equations for 
an asymmetric case £i ^ £ 2 . Clearly, the DPA model is a practical and broadly applicable 
tool for exploring the physics of a more complicated model. In the future we plan to use 
the DPA frame-work to derive, and to get physical insight in, an approximate dispersion 
relation for capillary waves on the interface, which takes into account the finite thickness 
and the structure of the interface. 
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